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ABSTRACT 

In  this  paper  we  carry  out  an  asymptotic  analysis  of  the  system  of 
differential  equations  describing  the  transient  behavior  of  a  p-n-junction 
device  (i.e.  a  diode).  We  determine  the  different  time-scales  present  in  the 
equations  and  investigate  which  of  them  actually  occur  in  physical 
situations.  We  derive  asymptotic  expansions  of  the  solution  and  perform  some 
numerical  experiments. 
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SIGNIFICANCE  AND  EXPLANATION 

The  present  paper  is  concerned  with  a  mixed  elliptic-parabolic  singularly 
perturbed  initial  boundary  value  problem  describing  a  p-n- junction  device. 

The  purpose  of  transient  modelling  is  the  determination  of  the  response  time 
of  a  particular  device  (as  opposed  to  the  stationary  model  where  the  emphasis 
lies  on  the  computation  of  voltage  current  characteristics).  We  carry  out  an 
asymptotic  analysis  using  a  singular  perturbation  parameter  which  is 
proportional  to  the  minimal  Debeye  length. 

Essentially  we  deal  with  the  involved  equations  on  2  different  time 
scales:  One  is  proportional  to  the  size  of  the  device.  The  other  one  is  the 
so  called  "dielectric  relaxation  time"  which  is  several  orders  of  magnitude 
smaller  and  essentially  a  property  of  the  material  and  the  doping 
concentration.  Thus  only  effects  occurring  on  the  "slow"  timescale  (which  is 
on  the  order  of  10  sec)  can  be  accelerated  by  miniaturization  of  the  device 
whereas  the  dielectric  relaxation  time  represents  a  lower  bound  to  this 
(technological)  approach.  This  fact  has  been  known  by  engineers  for  several 

years.  However,  this  analysis  gives  a  mathematical  proof  for  it.  Moreover, 

* 

it  is  precisely  determined  which  effects  are  due  to  which  time  scales.  We 
derive  asymptotic  expansions  for  the  solution  in  powers  of  the  perturbation 
parameter  and  give  an  existence  proof  for  the  solution  of  the  reduced 
problem.  We  determine  (and  prove)  the  decay  rate  of  the  fast  time  scale 
solution. 
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The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


AN  ASYMPTOTIC  ANALYSIS  OF  A  TRANSIENT  p-n  JUNCTION  MODEL 


Christian  Ringhofer 


0.  Introduction 

In  this  paper  we  present  an  asymptotic  analysis  of  a  singularly  perturbed  mixed 


system  of  parabolic  and  elliptic  equations  modelling  a  p-n-junction  device*  We  consider 
the  following  physical  situation:  A  semiconductor  (for  instance  silicone)  is  doped  with 
donor  atoms  (positive  ions)  in  the  left-hand-side  and  with  acceptor  atoms  (negative  ions) 
in  the  right-hand-side  of  the  device.  The  device  possesses  2  ohmic  contacts  to  which  a 

_3 

bias  is  applied.  The  devive  is  assumed  to  have  a  characteristic  length  21 (~  0.5*10  cm). 

The  physics  of  a  p-n  junction  are  explained  in  Sze  [1969],  Ashcroft  et  al.  [1976]  and  R. 

A.  smith  [1978].  The  equations  governing  the  potential  distribution  and  the  carrier  and 

current  densities  are  in  the  case  of  one  space  dimension  (see  Van  Roosbroeck  [1950]): 


♦  ”  ?  (n  -  p  -  C) 

xx  e  * 

n.  *  1  J  -  R 


Z  q  n 


J_  -  R 


Poisson's  equation 
continuity  equations 


for  electrons  and  holes 


Jn  -  <*(Dnnx  "  VV 
Jp  =  -qtDpP*  +  UpPV 


current  relations 


for  electrons  and  holes  . 


The  system  (0.0)-(0.4)  is  subject  to  the  boundary  conditions 
(0.5)  (a)  <M-i,t)  -  -^(t),  (b)  i|»U,t)  -  ^(t) 

(0.6)  n(-i,t)  *  p_i  +  C(-*>,  p(-i,t)  «  p_, 

(0.7)  n(i,t)  *  p,  +  C(£),  p(t,t)  ”  p, 


and  the  initial  conditions 


*(x,0)  »  ifijtx) 


n(x,0)  *  nj(x) 
p(x ,0 )  =  Pj(x)  . 
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The  meaning  of  the  dependent  variables  with  units  is  given  in  Table  1 : 


Table  1 

^  electric  potential  (V) 

electric  field  (V  cm-1) 
n  electron  density  (cm-3) 

p  hole  density  (cm-3) 

_2 

Jn  electron  current  density  (ft  cm  ) 

_2 

Jp  hole  current  density  (A  cm  ) 

The  parameters  q,  e,  un#  Up,  Dn,  Dp/  n^  UT  in  (0.0)  —  <0.4)  have  the  following  meaning 
and  approximate  value  at  T  »  300K  (room  temperature). 

Table  2 


<3 

elementary  charge 

10“19  As 

e 

permittivity  constant 

10~12  AsV- 3 cm" 

Un 

electron  mobility 

103  cm2V~1s“1 

U 

p 

hole  mobility 

103  cm2V~1s_1 

ni 

intrinsic  density 

1010  cm-3 

Dn 

electron  diffusion  constant 

25  cm2*”1 

D 

n 

hole  diffusion  constant 

25  cm2s_1 

P 


C(x)  in  (0.0)  is  a  given  function  of  space  and  models  the  doping  profile  (i.e.  the 
preconcentration  of  electrons  and  holes).  We  assume  that  C(x)  <  0  for  x  <  0  (in  the 
p-region)  and  C(x)  >  0  for  x  >  0  (in  the  n-region  holds).  Thus  we  locate  the  p-n 
junction  at  x  =  0.  Moreover  we  assume  that  C(x)  has  a  jump-discontinuity  at  x  »  0 
(abrupt  junction)  and  that  C(x)  is  an  odd  function  of  x  i.e. 

(0.9)  C(-x)  =  -C ( x)  for  x  f  0,  C(0)  f  0 
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holds.  (Although  the  oddness  of  C  does  not  appear  as  a  special  case  from  a  physical 
point  of  view  it  will  simplify  the  analysis  considerably).  Because  of  the  physical 
background  of  the  boundary  conditions  <0.5)— {0.7)  the  oddness  of  C(x)  results  in 
p_^  =  p^  +  C(l).  For  the  recombination  rate  R  in  (0.1),  (0.2)  we  take  the  so  called 
Shockley  Read  Hall  recombination  terms 


n*p  -  n. 


(0.10) 


Tn(n  +  nA  )  +  r  (p  +  nA) 


t  and  t  denote  the  average  lifetime  of  electrons  and  holes.  Realistic  numerical 
n  p 

values  for  t  and  r  are  t  »  t  »  10“®(s).  Other  recombination  rates  can  be  found 
n  p  n  p 

in  linger  et  al.  [1981]  and  Schiitz  et  al.  (1981).  Generally  it  can  be  said  that  the 
choice  of  the  recombination  rate  is  not  very  important  for  the  behavior  of  the  solution  of 
the  time  dependent  problem  (0.0)-(0.8)  as  long  as  we  do  not  consider  impact  ionization 
effects  (whereas  it  is  of  considerable  importance  in  the  steady  state  case  (see  Mock 
[1983])). 

For  the  rest  of  this  paper  we  assume  that  the  initial  functions  and  pj  in 

(0.8)  are  compatible  with  the  boundary  conditions  (0.5 )  —  <0.7)  and  the  (elliptic)  equation 
(0.0).  Thus  we  request  that 

(0.11)  »  f  (nx  -  Pl  -  C) 

XX 

(0.12)  p(±*)  -  P±1 

(0.13) 

(0.14)  ♦j(-l)  “  1  ( 0 ) ,  tjU)  -  11^(0) 

holds.  After  an  appropriate  scaling  we  carry  out  a  singular  perturbation  analysis  of 
problem  ( 0.0 )-( 0.8)  where  the  quantity 


n(±i )  •  p±1  +  C(±*  ) 


(0.15) 


X  -  (- 


E  U_ 


1/2 


1  q  max|C(x) | 
x 


(which  is  proportional  to  the  minimal  Debeye  length)  acts  as  a  perturbation  parameter.  In 
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this  analysis  we  will  rely  on  a  similar  asymptotic  analysis  for  the  steady  state  problem 

by  Markowich  and  Ringhofer  [1981]  and  Markowich  et  al.  [ 1982 J. 

It  turns  out  that  in  addition  to  the  internal  layer  behavior  already  present  in  the 

steady  state  case  (see  Markowich  and  Ringhofer  [1981]  there  exist  two  different  time 

scales  in  the  problem  (0.0)-(0.8).  One  time  scale  (the  "slow"  scale)  is  proportional  to 
2  —Q 

l  /u  and  is  of  the  order  of  magnitude  10  sec  for  modern  devices .  Here  u  ■  W  “  U 

n  p 

denotes  the  mobility.  The  other  ("fast")  time  scale  is  proportional  to  the  so  called 
dielectric  relaxation  time  (see  Sze  (1969])  which  is  several  orders  of  magnitude 
smaller.  Thus  the  asymptotic  solution  of  (0.0)-(0.8)  will  consist  of  4  different  parts. 
Two  parts  vary  on  the  slow  time  scale  -  one  of  them  valid  near  the  p-n  junction  (x  “  0) 
one  of  them  valid  away  from  x  ■  0.  The  two  other  parts  -  again  an  outer  solution  and  a 
p-n-junction  layer  solution  -  vary  on  the  fast  time  scale  and  are  only  valid  near  t  ■  0. 
From  the  work  of  Markowich  and  Ringhofer  [1982]  we  conclude  that  the  fast  time  scale  is 
not  present  if  'll-,  nt  and  pj  are  solutions  of  the  steady  state  problem 


(0. 16) 

XX 

.  a 

e 

(n  -  p  -  C) 

(0.17) 

(a) 

J  -  qR, 

X 

(b) 

Jn  “  ’'Vx  -  VV 

(0.18) 

(a) 

Jp  ■  -3R' 

Fx 

(b) 

Jp  -  -q<t>pPx  +  upp*x>  . 

Together  with  the  boundary  conditions  corresponding  to  (0.5)-(0.7)  in  the  steady  state 
case. 

This  paper  is  organized  as  it  follows:  In  Section  1  we  scale  (0.0)— (0.8)  and 
reformulate  the  problem  as  a  singular  perturbation  problem.  In  Section  2  we  derive  the 
asymptotic  expansion  for  the  slow  time  scale.  We  prove  the  existence  of  a  solution  of  the 
reduced  problem  and  give  necessary  and  sufficient  conditions  for  the  slow  time  scale 
expansion  to  satisfy  the  differential  equations  the  boundary  and  the  initial  conditions. 
These  conditions  are  satisfied  if  the  initial  data  i|i  ,  n ^  and  Pj  are  solutions  of  the 
steady  state  problem  ( 0. 16)-( 0. 18) .  In  Section  3  we  consider  the  case  when  these 
conditions  are  violated  and  derive  an  asymptotic  expansion  for  the  solution  of  (0.0)-(0.8) 
with  general  initial  data  (0.8)  which  just  satisfy  the  compatibility  conditions  (0.11)- 
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sss 


(0.14).  In  Section  4  we  present  some  asymptotic  solutions  which  were  obtained  by  solving 
the  reduced  problem  and  the  layer  equations  numerically  by  finite  difference  methods.  The 
proofs  of  some  results  stated  in  Sections  2  and  3  are  given  in  the  appendix. 


1.  Reformulation  as  Singular  Perturbation  Problem 

First  we  scale  (0.0)— (0.8)  appropriately  and  reformulate  it  as  a  singular 

perturbation  problems  We  scale  the  space  variable  x  by  the  characteristic  length  l 
2 

and  the  time  by  4  /D  where  D  «*  D_  =  D  denotes  the  diffusion  coefficient 

n  p 

(1.0)  *s=T'  ta=tV 


For  n,  p,  Jn,  J  we  use  the  following  scaling: 


TJ& 


't'  -  n  =  Cng,  p  =  Cp8  , 


(1.2)  J  =  J  ,  J  »  ^  J 

n  2  ns  p  2  ns 

(1.3)  C  j=  max|C(x) |  . 

x 

(Here  the  variables  with  the  subscript  s  denote  the  scaled  variables.) 

(1.4)  Remark :  For  the  rest  of  this  paper  we  assume  the  Einstein  relations 


D  D 
T  Wn  Up 


UT  (the  thermal  voltage)  is  only  a  function  of  the  temperature  (which  is  assumed  to  be 
constant ) . 

After  scaling  the  system  (0.0)— (0.8)  assumes  the  form 


,2  3  ,  „  „  ,  .  C(x) 

,2  ^s  "s  ~  ps  "  Cs'  cs(xs)  max | C ( x ) | 


<3)  3T  "s  *  JT  J"  “  V  (b)  Jn  “  3“  ST*  *. 


St  3X  “Pa  "s' 


n  3x  s  "s  3x 


'P.  3x  *8  *-8  3x  rs 

8  s  s 


n  p  -  6  n. 

_  _  s  s _  .  _  i 

Rs  _  n  +  p  +  26  '  6  - 

s  *s  C 


boundary  conditions: 


(1.10) 

V-1,t)  “  <t), 

V^V  "  *1  (t> 

s 

s 

(1.11) 

nst-1'ts)  ”  P-18  +  c(_1)' 

nS(1'tS)  “  p18  +  c<1) 

(1.12) 

Ps'^'V  “  P-18' 

Pst’/tg)  -  p,8 

initial  conditions 

(1.13) 

Vxs'°>  ■  ’•’l  (xs)'  ns(xs'0)  “ 
s 

nI8<xs*'  Ps{xs'0)  "  Plg^s'  • 

Here  i)>±1  , 

*!  ,  P±1  ,  "I  ,  Pl  denote  ip  , « 

S  1  8  8  8 

pj,  p±1»  nj,  pj  scaled  in  the  same  way  as 

<ti,  n  and 

form 

p.  The  (dimensionless)  singular  perturbation  parameter  X  in  (1.6)  is  of  the 

«0_ 

(1.14) 

x  -  (-JJ  . 

£  *  max |C(x) |  . 

qlC 

X 

A  is  proportional  to  the  minimal  Debeye  length  (see  See  [1969]).  From  hereon  we  will 
omit  the  subscript  8.  Since  C(x)  is  an  odd  function  we  can  employ  the  "Ansatz" 

(1.15)  ^(-x)  »  ->Mx),  n(x)  -  p(-x),  Jn(x)  -  Jp(-x) 

to  obtain  a  problem  posed  on  [0,1].  This  enables  us  to  treat  the  internal  layer  at  the 
p-n  junction  (see  Markowlch  and  Ringhofer  [1981])  as  a  boundary  layer.  In  order  to 
employ  this  "Ansatz"  we  have  to  assume  the  corresponding  relations  for  the  intial  data 
♦x*  nx  and  pj. 

(1.16)  Assumption:  The  initial  data  satisfy 

(1.17)  Pjl-x)  -  -^(x),  nj(x)  -  Pjl-x)  . 

After  employing  (1.15)  we  obtain  (1.6)“(1«13)  in  its  final  form: 

(  1.18)  X2^  -  n  -  p  -  C 

(1.19)  (a)  nt  -  Jn  -  R,  (b)  J„  -  nx  -  ni|>x 


(1.20)  (a)  pt  -  -Jp^  -  R,  (b)  Jp  -  -px  -  p*x 


(1.21) 


boundary  conditions: 


n  •  p  -  5 
n  +  p  +  2i 


(  1.22) 


(1.23) 


initial  conditions: 


(1.24) 


♦  (1«t>  =  tl^tt),  n(1,t)  “  p1  +  C(  1),  p(1,t)  -  p, 
<H0,t)  -  0,  n(0,t)  -  p(  0, t )  ,  Jn(0,t)  -  Jp(0,t) 

<Hx,0)  -  ^jfx) ,  n(  x,  0 )  -  nx(x),  p(x,0)  -  Pj(x)  . 


(1.25)  Notation:  For  the  rest  of  this  paper  we  denote  the  solution  vector 


<i|i,n,p,J  ,J  )  by 
n  p 

(1.26) 


w  -  (i(i,n,p,J  ,J  ) 
n  p 


The  boundary  conditions  (1.23)  emerge  from  the  "Ansatz"  (1.15).  The  compatibility 


conditions  now  read 


(1.27) 


(1.28) 


(1.29) 


X  Vj  “  ni  "  Pi  "  c 

XX 

*j(1)  -  <!>,(*),  nj(1)  -  b,  Pj(1)  -  a 
<1^(0)  -  0,  nj(0)  -  Pj(0),  nT  (0)  -  -Pj  (0) 


Again  (1.29)  is  obtained  from  (1.17)  and  guarantees  (1.23)  to  hold  at  t  *  JB.  In  practice 
♦  nj  and  px  will  be  the  (scaled)  solution  of  the  steady  state  problem  ( 0. 16 )  —  (  0. 18) . 
From  Markowich  and  Ringhofer  [1981]  we  know  that  this  solution  has  an  asymptotic  expansion 
in  powers  of  X  consisting  of  an  outer  solution  and  layer  terms  varying  in  the  variable 
5  “  x/X.  From  here  on  we  will  assume  that  ’t'j,  nx  and  px  have  an  asymptotic  expansion 
of  the  same  form: 


(1.30) 


Wj(x)  «  w  (x,X)  ~  l  [w3(x)  +  wr(f)lx:1 
j-0 


WI  “  (l|'i'ni,pi> 


where  the  w^  satisfy 


(1.31) 


2 ' 

*Wj(£)l  <  C^e  for  some  constants 


\y:  * . 

k.V.V.V, 
k  V  V  V 

4  "r*.  V  1 

h ^  •  ■  a** 

A  sVwV. 


k! 


Moreover,  we  will  assume  that  the  initial  data  are  such  that  the  current  densities  Jn 
and  Jp  stay  uniformly  bounded  at  t  “  0  as  X  tends  to  zero: 


(  1.32) 


(  1.33) 


x.0,1]  n 


J  (x,X)|,  max  |j  (x,X)|  <  const,  as  X  ♦  0 


xe[0,1] 


J  _  nT  -  n*  , 


“Pi,  -  PI*I 


2.  Slow  Time  Scale  Expansions 

In  this  section  we  derive  our  asymptotic  expansion  to  the  solution  w  of  (1.18)- 
(1.24)  of  the  form 


w(x,t,X)  ~  £  (w:'(x,t)  +  w-*(^,t)]X^ 

j  =  0 


Thus  we  obtain  an  approximation  of  the  solution  on  the  "slow"  t-time  scale.  The  expansion 

(2.1)  will  satisfy  the  differential  equations  ( 1. 1 8 ) — ( 1.20)  and  the  boundary  conditions 

( 1.22)— (1.23).  However,  it  will  not  satisfy  the  initial  conditions  (1.24)  for  arbitrary 
initial  data  i|>  ,  nt  and  pj.  We  derive  a  necessary  and  sufficient  condition  on 
(ifr^n  ,p  )  such  that  the  expansion  (2.1)  satisfies  (1.24)  as  well.  As  it  turns  out 
these  conditions  are  satisfied  if  n^  and  p^  in  (1.24)  are  the  solution  of  the 

steady  state  problem 

(2.2)  X%xx  -  n  -  p  -  C 

(2.3)  -R  '  Jnx'  Jn  “  nx  ’  n*x 

(2.4)  R  “  /  J  »  -P„  -  P+ 

*X  P  * 

together  with  the  boundary  conditions  ( 1.22)-( 1.23).  Since  this  will  be  the  case  in  all 
"physically  meaningful"  situations  (2.1)  actually  represents  an  approximation  to  the 
solution  of  ( 1. 18)-( 1.24).  Furthermore  we  give  an  existence  result  for  the  zero'th  order 
term  w°  in  (2.1). 

We  start  with  the  equations  determining  w° :  Setting  \  equal  to  zero  in  (1.18)- 
(1.24)  gives 

„  ~0  ~0 

(2.5)  0  «  n  -  p  -  C 


~0  ~0  ~0 
(a)  n  =  J  +  R  , 
t  n 

x 


~0  ~0  -0-0 

(b)  J  *  n  -  n  if 
n  x  x 


~o  ~o  ~o  <LV  ~o  -0  -o-o 

(a)  pt  =  -J  +  R  ,  (b)  J  =  -px  -  p 

Px  P 


-0  -0  ,2 
n  •  p  -  6 


n°  +  p°  +  25 


Thus  Poisson's  equation  (1.18)  degenerates  to  the  condition  of  vanishing  space  charge  for 


X  ♦  9. 


Subtracting  (2.7)  from  (2.6)  and  using  (2.5)  we  can  rewrite  (2.6)-(2.7)  into 


t(2P°  +  C)*x]x  *  Cxx 


(2.10) 


~o  ~0  ~o  ~o  ~o  ~o~o 

(a)  p  -  -J  +  R  ,  (b)  J  -  -p  -  p  i{i 

t  p  p  x  x 

X 


(2.11) 


~0  ~0 

~0  p  (o  +  C 
R  *  r  - 

2p°  +  C 


(2.9)  is  an  elliptic  equation  expressing  the  conservation  of  the  total  current  Jn  ♦  Jp 


(in  the  absence  of  space  charge)  coupled  to  the  parabolic  conservation  law  ( 2. 10) (a) ,(b) . 


~0  ~0 

n  and  J  are  then  given  a  posteriori  by  virtue  of  (2.5): 
n 

~0  ~0 

(2.12)  n  -  p  +  C 


(2.13) 


~0  ~0 


”  Px  +  °X  "  (P  + 


Since  (2.9)-(2.10)  is  a  system  involving  two  second  order  differential  operators  in  space 


we  can  not  expect  its  solution  to  satisfy  all  six  boundary  conditions  ( 1.22)— ( 1.23) .  At 


the  contact  (x  -  1),  however,  the  boundary  conditions  have  been  derived  from  the 


assumption  of  vanishing  space  charge  (n  -  p  -  C  *  0).  Therefore  the  boundary  conditions 


(  1. 22)-(  1 . 23)  are  consistent  with  (2.5):  If  we  impose  the  boundary  conditions 


(2.14) 


*°(1,t)  -  ^(t),  p°(1,t)  -  pn 


~0  ~o 

n  :«  p  +  C  satisfies  automatically  the  third  boundary  condition 


(2.15) 


n  ( 1 ,t)  ■  p1  +  C( 1 )  . 


Therefore  no  boundary  layer  term  is  needed  at  x  ■  1  (the  contact).  At  x  -  0  (the 


~0  ~0 

junction)  the  condition  n(0,t)  »  p(0,t)  contradicts  n  -  p  -  C  ■  0.  Therefore  the 


outer  solution  w  has  to  be  supplemented  by  a  boundary  layer  term  there.  (In  other 


words  the  approximation  of  vanishing  space  charge  can  not  be  valid  near  the  p-n 


junction.)  Thus  we  supplement  w  (x,t)  by  the  layer  term 


(2.16) 


x  “0  “0  '0  *0  10  10. 


w  (S,t),  5  »  T  ,  w  «  0|»  ,n  ,p  ,J  ,J  ) 
a  n  p 


where  we  request  that  Iw  8  decays  exponentially  as  £  ♦  Inserting 


*\  •*,  r .  * « 

v  v"v  V 


.  *  .  '  a  *  .  ’  -  ' 

1  .  e  *  »  *  .X»  ‘  .  ‘ 

*.s 


. 


w  (X£,t)  +  w  (£,t)  into  ( 1 . 1 8  >  —  ( 1. 20)  and  X  *  0  gives  (using  (2.9)-(2. 13) ): 


(2.17) 


(2.18) 


(a)  0  =  J  , 

nc 


“0  “0  *0 
*  n  -  P 

*0  ~0  *0  *0 
(b)  0  -  nj.  -  (n  (0,t)  +  n  )\|)^ 


(2.19) 


(a)  0  -  -J  ,  (b)  0  -  -p^  -  (p  ( 0,t)  +  p  )i|>£ 


Because  w  has  to  vanish  at  £  ”  "  we  can  integrate  ( 2. 18)— ( 2. 19)  and  obtain 


(2.20) 


(2.21) 


0  *0 

J  =  0,  J  =  0 
n  P 

n°(C,t)  =*  n°(0,t)  [exp(\|i°(?»t) )  -  1] 


(2.22) 


p°(5,t)  -  p°(0,t)  [exp(-i(i0(C,t))  -  1] 


Inserting  w  +  w  into  the  boundary  conditions  (1.23)  and  X  ♦  f  gives 


(2.24) 


*0  ~0 
♦  (0,t)  =  -i|>  (0,t) 


(2.25) 


n°(0,t)exp(i(i0(0,t))  -  p°(0,t)exp(-i|'°(0,t)) 


(2.26) 


~0  ~0 

J  (0,t)  =  J  (0,t)  . 
n  p 


From  ( 2 . 24 ) - ( 2 . 26)  we  obtain  two  boundary  conditions  for  the  system  (2.9)-(2.10)  at 


x  ■  0.  Thus  we  have  the  following 


(2.27)  Theorem!  The  zero'th  order  term  w  of  the  outer  solution  satisfies  the  reduced 


problem: 


(2.28) 


((29°  +  C)*°]  -  C 

r  x  x  xx 


(2.29) 


~0  ~0  ~0 
(a)  Pt  -  -J  +  R  , 

px 


..  .  ~0  ~0  ~0~0 

(b)  J  «  -p  -  p  ♦ 
P  X  X 


(2.30) 


(p°(0,t)  +  C(0)]exp[-9°(0,t)]  »  p°(  0,t)expU>°(  0,t)] 


(2.31) 


(2.32) 


~0  ~0 

J  (0,t)  -  J  (0,t) 
n  p 

*°<1,t>  -  ♦1(t),  p°(1,t)  -  p, 


(2.33) 


~o  ~0  .  _ 
n  “  p  +  C, 


~0  ~0  ~0  ~0 
J  -  p  +  C  -  p  +  C)<i 
n  x  x  x 


(2.34) 


p°(x,0)  =  f ( x ) 


v'- •.*%-;*.• V\\V\  ■.■'.•VW./%V’>WAv  V.v.v.v. v.v. 


-  ■  VV  •  V 


(2.35)  Remark;  We  have  not  determined  the  initial  function  f(x)  yet.  The  appropriate 
choice  of  f(x)  will  be  discussed  at  the  end  of  this  section  and  in  Section  3.  In  any 
case  we  can  prove  an  existence  result  for  small  f  up  to  a  certain  value  of  the  bias 

(2.36)  Theorem;  There  exist  positive  constants  a.  and  K  such  that  if 


(2.37) 


a  ■  max|f(x)|  <  and  \|>  >  tn  — 


holds  then  ( 2. 28)-( 2. 34)  has  a  classical  solution  which  is  continuously  differentiable 
with  respect  to  x  and  t  for  x  e  10,11  and  all  t  less  than  some  arbitrary  but 

finite  time  T.  The  proof  of  Theorem  (2.36)  is  deferred  to  the  Appendix. 

~0  *0 
Once  w  is  obtained  as  the  solution  of  problem  (2. 28)- (2. 34)  w  can  be  obtained 

by  solving  (2. 17)— <2. 19)  together  with  the  appropriate  boundary  conditions: 

“0  *0  *0  *0  ‘0  *0 

(2.38)  Theorem:  The  layer  term  w  -  (i(i  ,n  ,p  ,J  ,J  )  in  (2.1)  satisfies 
-  n  p 


(2.39) 


(2.40) 


*0  ~0  ill0  ~  -i(»0 

iji^j.  -  nu(0,t)[ev  -  1]  -  p(0,t)[e  v  -  1] 

*0  ~0  *0 

(0,t)  -  -1(1  (0,t),  III  (»,t)  -  0 


(2.41)  n°(E,t)  -£°(0,t)te*  (C,t>  -  1),  p0(E,t)  -  p\o,t)  [e-*  (5,t>  -  1] 


(2.42) 


jV,t)  -  J°(5.t>  »  0 

n  p 


We  now  turn  to  the  question  of  an  appropriate  choice  for  the  initial  function  f(x) 
in  (2.34).  As  mentioned  earlier  we  can  not  satisfy  the  initial  conditions  (1.24)  with  our 
expansion  for  arbitrary  initial  functions  which  just  satisfy  the  compatibility 

condition  (  1.27)-(  1.29) .  From  the  previous  asymptotic  analysis  we  conclude 

(2.43)  Lemma:  If  the  solution  w(x,t,A)  of  ( 1 . 18)-( 1 . 24)  possesses  an  asymptotic 


expansion  of  the  form  (2.17)  then  the  zero'th  order  term  Wj(x)  in  (1.30)  has  to  satisfy 


(2.44) 


t(2p°  +  C>*°  ] 

x  x 


~0  *0  . 

If  (2.44),  is  satisfied  then  the  zero  order  term  w  (x,t)  +  w  (E,t)  in  (2.27),  (2.38) 
satisfies  the  differential  equations  the  boundary  and  the  initial  conditions  up  to  terms 
of  order  0(A)  if  the  Initial  function  f(x)  in  (2.34)  is  chosen  as  p°(x). 


Proof :  Since  (2.28)  is  time  independent  it  obviously  has  to  hold  for  t  “  0  if 

w(x,t,X)  has  an  asymptotic  expansion  of  the  form  (2.1).  Therefore  (2.44)  is  a  necessary 

~0  *0 

condition.  On  the  other  hand  w  (x,t)  +  w  (C.t)  satisfy  the  differential  equations  and 
the  boundary  conditions  up  to  terms  of  order  0(1).  Hie  compatibility  condition  (2.46) 
impl ies 


(2.47) 


0 


-  C 


(2.48) 


P 


0 

I 


(2.49) 


'0 

1^(0) 


-J°(0) 


If  we  choose  f(x)  in  (2.34)  to  be  p^(x)  then  (2.33)  and  (2.27)  imply 

~0  ~0  ~0 
n  (x,0)  =  n^(x).  In  this  case  <|)  (x,0)  satisfies 


(2.50) 


( (  2p  J  ( x )  +  C(x)  )*°(x,0))x  -  Cxx(x) 


(2.51) 


lp°  (0 )  +  C(0)]e"*  (0'0)  -  Pj(0)e^  <0'0),  *°(1,0)  -  ♦1  (0)  . 


From  (1.32)  we  conclude 

(2.52)  J°  -  nT  -  <n_ ( 0 )  +  nT>*  -  0 

"i  h  1  1  H 


since  n^  decays  exponentially  for  £  ♦  • 

(2.53)  n^O  -  nI(0)[e  1  -  1] 
and  similarly 

"*I 

(2.54)  Px  (C )  -  pt(0)  [e  -  1] 

has  to  hold.  Inserting  +  '  nj  +  nj  >  Pj  +  Pj  into  (1.29)  and  X  ♦  JJ  gives  together 


with  (2.47) 


b* 


s. 

I 

f. 


f. 

r 

f 

HU 

r 


I 


p5 


(2.55) 


-’t'j(O)  _Q  tj(O) 

-  V 


Thus  and  i(>(x,0)  satisfy  the  same  differential  equation  (2.50)  and  the  same  boundary 

conditions.  Since  this  problem  has  a  unique  solution  (which  follows  from  a  simple  maximum 

principle)  >Ji°(x)  =  ^°(x,0)  has  to  hold.  Inserting  (2.53)  and  (2.54)  into  (2.48)  gives 

*0  '0 

the  same  boundary  value  problem  for  ^  as  <(i  (5,0)  has  to  satisfy  namely 


(2.56) 


(2.57) 


a  d;  — 

A0  ~0,B4,  yI  ,,  ~0...  .  yI 

'll  »  n  (0)  (e  -  1]  -  p  (0)[e 

*55  I  I 


4»°( o )  -  -JjfO),  if^f*)  “  0  . 


1] 


A0  *0 

Since  (2.56)-(2.57)  again  has  a  unique  solution  4^(5)  “  ♦  (5,0)  and  therefore 
n°(5)  ■  n°(5,0)  and  p°(5)  ”  p°(5,0)  has  to  hold. 

(2.58)  Remark!  So  if  the  initial  functions  *  ,  nx  and  px  satisfy  (2.44)  in  addition 
to  (2.58)  we  have  constructed  a  zero  order  asymptotic  solution  to  problem  (  1. 18)  —  ( 1.24) . 
Higher  order  expansions  can  now  be  obtained  in  a  straightforward  manner.  If  the 
corresponding  assumptions  are  made  on  the  higher  order  terms  of  tji^,  nx  and  px  these 
higher  order  expansions  satisfy  ( 1. 18 )  — <  1 . 24)  to  any  arbitrary  order  in  X. 

We  now  conclude  this  section  with  the  following 

(2.59)  Theorem;  If  the  initial  functions  i|>  ,  nx  and  px  are  the  solution  of  the  steady 
state  problem 


(2.60) 

X\ 

-  nx  -  Pj  -  C 

XX 

(2.61) 

(a)  J  -  R, 

Ax 

(b)  J  -  nx  -  n^j 

lx  x 

(2.62) 

(a)  J  »  -R, 

Plx 

(b)  J  -  -Pj  -  Pj* 

1  *  X 

(2.63) 

^(0)  -  0  , 

♦jtl)  -  (0) 

nI<0>  “  Pl(0)  • 

nx( 1 )  -  n, 

t. 

C. 

?. 
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(2.64) 


(0)  -  (0)  , 


Pl(D  -  Pi 


then  the  conditions  (2.44)  is  satisfied. 

For  the  proof  of  this  statement  the  reader  is  referred  to  Markovich  and  Ringhofer 
(1982)  where  a  singular  perturbation  analysis  of  the  problem  (  1. 18)-(  1.24)  can  be  found. 


V  \ 


In  this  section  we  investigate  the  asymptotic  behavior  of  the  solution  of  (1.18)- 
(1.24)  if  the  conditions  ( 2.44)-( 2. 45)  on  the  initial  functions  are  violated.  It  turns 
out  that  in  this  case  there  is  an  additional  "fast"  time  scale  present  in  the  solution. 
Although  this  time  scale  will  not  be  present  in  the  physically  meaningful  situations,  as 
we  saw  in  the  preceding  section  (since  in  these  cases  the  initial  functions  are  solutions 
of  the  steady  state  problem  (2. 17)-( 2. 19) ) ,  it  is  of  some  interest  to  investigate  these 
"fast"  solutions. 

Any  perturbation  of  the  solution  will  propagate  on  this  time  scale.  Therefore  the 
understanding  of  this  fast  time  scale  behavior  is  important  for  the  development  of  any 
numerical  method  for  the  solution  of  ( 1. 18)— { 1.24)  (see  Ringhofer  [1983]).  Thus  we  derive 
our  asymptotic  solution  of  ( 1. 18)-(  1.24)  for  arbitrary  initial  functions  nj  and  Pj 

which  just  satisfy  {  1 . 27 )  —  (  1.29)  (and  not  necessarily  ( 2. 44)-( 2.45) ) . 

First  we  determine  the  possible  fast  time  scale.  For  this  purpose  it  is  convenient 
to  reformulate  problem  ( 1 . 1 8 ) — ( 1 . 24 ) ;  Differentiating  Poisson’s  equation  with  respect  to 
time  yields 

(3.D  ^  "  nt  -  Pt  • 

2 

Inserting  for  nt  and  pt  from  ( 1. 19)-( 1.20)  and  replacing  n  by  p  +  C  ♦  X  f  gives 

(3.2)  (a)  -  Jx,  (b)  J<-  ♦  V  -  cx  ♦  *\xx  ~  (2P  ♦  c  ♦  X 

(3.3)  (a)  pt  -  Jp^  -  R,  (b)  Jp  -  -Pjt  -  p*x  . 

n  and  Jn  are  then  given  a  posteriori  by 

(3.4)  n  •  p  +  C  ♦  *2'l' 

xx 

(3.5)  Jn  “  J  ”  Jp  “  nx  ”  nl^x  ’ 

The  boundary  conditions  (  1 . 22 )-(  1. 23)  become 

(3.6)  *(0,t)  -  0,  X2*xx<0,t)  -  C( 0 ) ,  J(0,t>  -  2Jp( 0,t) 

(3.7)  *(1,t)  -  <!>,( t ) ,  X24»xxC  )  -  0,  p(1,t>  -  p,  . 
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Obviously  (3.2),  (3.3),  (3.6),  (3.7)  are  equivalent  to  ( 1- 18)  —  ( 1.24)  (provided  the 

solutions  are  sufficiently  smooth) . 

We  now  try  to  determine  a  possible  fast  time  scale: 

(3.8)  t  -  —  .  a  >  0  . 

X° 


Away  from  the  p-n  junction  at  x  -  0  where  we  expect  the  derivatives  w.r.t.x  to  be 
bounded  (3.2)(a)(b)  assumes  the  form 

(3.9)  x2*xxt  ’  Cxx  “  (2p  +  ^xx  "  (2p  +  C)x’*'x  +  0<x2)  ' 

If  2p  +  C  is  in  some  sense  "almost  constant"  away  from  x  -  0  (which  is  the  case  for  a 

constant  doping  profile  and  reverse  bias)  (3.9)  behaves  like  an  ordinary  differential 

2 

equation  for  and  therefore  exhibits  the  fast  time  scale  t/X  .  This  heuristic 

argument  can  be  made  more  rigorous:  We  supplement  the  asymptotic  expansion  (2.1)  by  an 
additional  term  varying  on  the  fast  time  scale  x  *  t/X°  where  o  >  0  has  to  be 
determined  yet: 


w(x,t,X)  ~  l  (w3(x,t)  +  wj(£,t)  +  *3(x,t)1 
j-0 


(3.10) 


5  -  |  ,  x  -  tX-01  . 

P  -  (♦3,u3,;3,3j,3fJ)T  . 

(w3  and  w3  are  the  terms  of  the  expansion  (2.1)  derived  in  the  previous  chapter.)  We 
assume  that  z3  decays  exponentially  as  T  +  ». 


(3.11) 


Then  we  have  the  following: 


-C2X 

Izj(x,t)I  <  C^e 


(3.12)  Lemma:  The  only  value  of  a  for  which  there  exists  a  non-trivial  solution  z  of 
the  corresponding  equations  is  a  -  2. 

Hie  proof  of  the  lemma  is  trivial  but  lengthy  and  is  therefor  deferred  to  the 
appendix. 


’  V  V  v  %  \  "  V  V*  ■•»•••  ,«•••*  •  *  i  fa*  *  «  «  .  •>  .  V.  «  .  r.  /  «.  » 


a 

.*  V 


(3.13)  Remark!  The  time  scale  t  «  t/X  has  the  following  physical  meaningi  In  unsealed 

form  it  is  of  the  size  eu  (DqC  )  1  sec  which  is  proportional  to  the  dielectric 

T  max 

relaxation  time  (see  Sze  [1969]). 

~0 

The  equations  for  z  are  now  obtained  from  (3.2).  Inserting  the  expansion  (3.10) 

2 

into  (3.2)— (3.9),  fixing  x  and  T  -  t/X  and  X  ♦  0  givesi 


~0  ~0  ~0  ~0  -^(W0 

(3.14)  *xxT  -  ~  [  ( 2p  +  2v  +  C)*x  -  2v  I|»xlx 

~o  ~0  ~0  ~0  ~0  ~0  ~0~0 

(3.15)  (a)  vT  «  0,  (b)  »  -vx  -  (p  +  v  )$x  -  v 

~0  -0 

(3.16)  u  ■  v 

~0  ~0  ~0  ~0  ~0  ~0~0 

(3.17)  J  -  u  -(n  +  u  )♦  -  ui(i  . 

u  x  x  x 

~0  ~0  ~0  ~0  ~0 
Here  4>  ,  n  and  p  have  to  be  taken  at  (x,0).  From  (3.15),  we  conclude  u  "  v 

~0 

So  we  obtain  4  from  the  equation 


(3.18) 


*xxT  -  -[(2p°(x,0)  +  C(x )  )♦  x  • 


~ 0  -0  ~0  's*0 

u  ,v  ,Jy,Jv  are  then  given  a  posteriori  by 

-0  -0 

(3.19)  u  -  v  -  0 

(3.20)  (a)  J°(x,r)  -  -p°(x,0)*x(x,t) ,  (b)  J°(x,t)  -  -n°(x,0  >$x<x,T )  . 

Thus  there  is  (in  the  zero  order  approximation)  no  fast  time  scale  correction  to  the 

~o  ~o 

carrier  densities  n  and  p  (u  «  v  •  0).  This  is  reasonable  since  we  expect  the 

carrier  flux  to  be  much  slower  than  the  dielectric  relaxation  time.  Correspondingly  the 

fast  time  scale  correction  of  the  currents  only  involves  the  drift  current.  To  determine 

a  solution  of  (3.18)  it  is  necessary  to  impose  two  boundary  conditions  at  x  ■  J 

and  x  -  1  and  an  initial  condition  at  T  -  0.  Since  we  have  already  satisfied  all 

~0  *0 

boundary  conditions  with  the  slow  time  scale  solution  w  +  w  we  obtain 


(3.21) 


*0(0,T)«  *°(1,T)  -  t 


We  can  prove  the  following  result  about  the  decay  of  $  i 
(3.22)  Theorem;  The  solution  of  the  problem 


(3.23) 


(3.24) 


♦  °(0,T)  -  ft,  $0(1,T)  -  0,  $°(x,0)  -  g(x> 


satisfies 


%'-/./VVV VV\JV-.  V  '  •  v.v.v.v  .-  *...vVv 


(3.25)  |$°(x,T)|  <  Ce 

for  some  positive  constant  C  and  any  a  satisfying 

~0 

(3.26)  2p  (x,0>  +  C(x)  >  o  . 

The  proof  is  deferred  to  the  appendix. 

We  now  investigate  the  choice  of  proper  initial  functions  for  the  problems  (2.28)- 
(2.34),  ( 3. 23 ) — ( 3.24)  in  order  to  construct  an  approximation  to  the  solution  of  (1.18)- 
(1.24)  for  general  initial  data  ip  ,  nx  and  px  which  just  satisfy  the  conditions 

~0  ~o 

( 1. 27 )  — ( 1. 32) .  Since  the  zero-order-fast-time-scale  corrections  u  ,v  vanish  (see 

(3.19))  we  have  to  choose  p°(x)  as  an  initial  function  f  for  problem  (2.28)  in  any 

case.  This  leads  to  the  following  procedure  to  obtain  a  zero  order  approximation  to  the 

solution  of  ( 1. 1 S ) — ( 1.24)  for  arbitrary  Initial  data: 

~0  -0  ~0 

Step  1:  Solve  problem  (2. 28)-( 2. 34)  with  f(x)(«  p  (x,0)>  •  Pj(x).  Obtain  <»  (x,t)  and 

n°(x,t).  n°(x,0)  -  nj(x)  is  then  satisfied  automatically  because  of  (1.27). 

A  o  *0 

Step  2:  Solve  the  slow  time-scale-boundary-layer  problem  (2. 39)- (2. 40) .  Obtain  ♦  ,  n 

and  p°.  Because  of  conditions  ( 1.27)-(  1.33)  ♦°C C » 0 )  “  ♦^(5),  n°(£,0)  «  n!j!(fi)» 

p°(5#0)  -  p°<£)  holds. 

If  *  ,  nx  and  px  are  the  solution  of  the  steady  state  problem  we  are  done.  If  not  we 
have  to  perform 

Step  3:  Obtain  by  solving  ( 3.  2 3)  —  C 3 .  24)  with  g(x)  ■>  'pj(x)  -  <i°(x,0). 

Now  the  solutions  obtained  by  Steps  1-3  satisfy  the  differential  equations  together  with 
the  boundary  and  initial  conditions  up  to  terms  of  order  O(X).  Higher  order 
approximations  can  now  be  obtained  in  a  straightforward  manner. 
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4.  Numerical  Examples 

In  this  section  we  present  some  numerical  results  for  a  simple  testproblem:  We  chose 
a  diode  of  length  1u  and  a  doping  profile  C  which  is  constantly  lO^cm-^  on  the  n-side 
(and  because  of  symmetry  -10  cm  in  the  p-region) .  This  corresponds  to  X  «  0.05.  We 
varied  the  applied  bias  from  0.5V  reverse  bias  to  0.5V  forward  bias.  Figures  1-7  show 
the  solution  for  x  €  [0,1]  where  the  p-n- junction  is  located  at  x  »  0  and  the  contact 
at  x  “  1.  To  obtain  a  solution  on  [-1,1]  the  graphs  in  Figures  1-7  would  have  to  be 
continued  according  to  (1.15).  For  all  calculations  backward  finite  differences  have  been 
used  in  time.  In  the  x-direction  we  used  the  exponentially  fitted  finite  difference 
method  introduced  by  Schar fetter  and  Cummel  (see  Scharfetter  and  Gummel  [1969]  and 
Markowich  et  al.  [1983b])  with  an  exponentially  graded  mesh  near  x  *>  0  to  resolve  the 
junction  layer.  Figures  1  and  2  show  the  "full"  solutions  of  ( 1. 18)-( 1.24)  for  p  and 
<|>  respectively.  As  initial  data  the  steady  state  solution  for  0.5V  reverse  bias  has 
been  used.  Thus  no  fast  time  scale  solutions  occur  near  t  ■  0.  Figures  3  and  4  show  the 
corresponding  solutions  of  the  reduced  problem  ( 2.28)-( 2. 34) .  In  Figures  2  and  4 
log1Q(p)  instead  of  p  is  plotted.  To  examine  the  fast  time  scale  behavior  we  held  the 
bias  constant  at  )JV  (equilibrium)  and  perturbed  the  initial  conditions  (thus  violating 
(2.44)).  The  dotted  line  in  Figure  5  is  the  equilibrium  solution  for  4>.  The  other  line 

is  the  initial  condition  for  <*.  Figure  6  shows  the  solution  for  of  the  full  problem 

t  ~o 

( 1. 18 >  —  { 1.24)  on  the  fast  time  scale  — .  Figure  7  shows  the  initial  layer  term  $  (x,x) 

XZ 

in  (3.10). 
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FIGURE  2 


FIGURE  3 
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5.  Appendix 

(5.1)  Proof  of  Theorem  (2.36);  Let  o  be  given  by 

(5.2) 

He  employ  the  following  transformation: 


a  «  max|f(x) |  . 
x 


p(x,t)  -  ap( x,t) ,  Jp(x,t)  -  ajp(x,t) 


(5.3) 


*<x,t)  -  a$(x,t)  +  i^t)  +  tn(|^-) 

(5.4)  ( 2.  28)-( 2. 34)  then  becomes 


(5.5) 


(5.6) 


t(2ap  +  C) 4*^  +  2p(tn  C)xlx  -  0 
(a)  pt  -  -3p  -  5,  (b)  Zfp  -  -px  -  p[(tn  C)x  +  a*x] 


p(ap  +  c) 


2ap  +  C  +  2w 


(5.7) 


(5.8) 


(5.9) 


p(x,0)  »  f(x>,  f(x)  :• 


f(x) 


(a)  p< 0, t)  -  [ap(0,t)  +  C(0)]e‘2a’1’<0't>  •  p(a,t) 

(b)  2p  ( 0 ,  t )  -  C(0)5  (0,t> 

*  X 

pO.t)  “  P1  -  1(1),  5(1, t)  -  JR  . 


p(a,t)  is  dependent  on  the  size  of  the  initial  data  (a)  and  the  applied  bias  (i|ii)) 


(5.10) 


cm  “"V**  1 

o(a.t)  -  •  e  •  r 


«  a0)  . 


C(0)  -  a 

If  a  is  small  enough  and  is  such  that  p  is  bounded  (5.5)-(5.9)  is  up  to  a  small 

perturbation  linear  in  5  and  p.  Thus  it  can  be  solved  by  Picard  iteration  if  the 
appropriate  spaces  are  chosen,  we  define 

(5.11)  Y  { ( u , v )  |  u,v  ec2[0,1]  *  (0,T) ,  0(1, t)  -  0  ,  C(0)Uv(0,t)  -  2v  <0,t)} 


Z  (C[0,1J  x  [0,T]  )2  x  C2[0,T]  *  C2[0,1J 


and  the  operators 
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q.  *  *  <  .  *  .  . 


A(!|»,p)  «  (((OM  +  2p(fn  C)x)x,  Pt  -  Pxx  -  tpUn  C)x)x  ,  p(0,t),  p(1,t),  p(x,0)) 

. . .  .  ......  _  .  -2aiM0,t> . . 


G(*,p)  -  <-2a(pijMx,  P^x> x»  P(®»t)e 


(ap(0,t>  +  C(0)),  0,0) 


Problem  ( 5.5)— ( 5.9)  can  then  be  written  In  the  form 


A(*,p>  -  G^p)  -  ( 0, 0, 0 ,? ( 1 ) ,? ( x) )  . 


A  Is  a  linear  operator  from  Y  onto  Z  with  a  bounded  Inverse  A  •  (This  can  Bimply 


be  shown  by  using  Fourier  expansions.)  The  Frechet  derivative  G*  of  G  satisfies 


(5. 12) 


IG*(p,K)l  <  a((2C(0)  +  a|p(0) | )p(a,t)e  V  ] 


Thus  it  can  he  shown  by  the  usual  Picard  iteration  argument  that  the  sequence  defined  by 


(5.13) 


A<Vl'W  “  G(W  + 


converges  to  a  solution  (p,<l<)  of  {  5.5)  —  ( 5.9) . 


(5.14)  Proof  of  Lemma  (3.12)s  Inserting  the  expansion  (3.10)  into  (3.2)  (which  are 


exponentially  small  away  from  x  -  0)  gives 


(5.15) 


2-a~0  2~  ~0  2~0  ~0  .  2~0  ~0 

X  $  ■*  X  $  -  ( ( 2p  +  C  +  X  +  2U  +  X  $  )♦] 

XXT  xxxx  *  XX  XX  X  1 


~0  2~0  ~0 

-  [  ( 2U  +  X>  )<|i  ]  ?  h.o.t. 

XX  X  X 


(5.16) 


— 0  a  ~0  ~0  —  0  ~0  ~0~0  or- 0 

uT  -  X  tux  ♦  (p  +  u  >*x  ♦  0  Vx)x  +  X  R 


Here  if»°  and  p°  have  to  be  evaluated  at  (x,tX°).  (5.16)  gives 


0°(x,T)  -  0 


for  X  *  0.  If  a  >  2  holds  then  (5.15)  becomes 


(5.17) 


♦°  -0 
XXT 


for  X  0.  Together  with  the  boundary  conditions  (3.21)  this  gives 


(5.18) 


-0 

*  (x,T )  »  0 


But  in  this  case  z  =  0  follows  immediately.  If  a  <  2  holds  then  (5.15)  gives  for 


(5. 19) 


~0  ~0 

■U2p  +  C)*xJx  -  0 


Again  using  the  boundary  conditions  (3.21)  gives  only  the  trivial  solution  for  ♦  (and 


therefore  also  for  z  ).  Thus  a  -  2  has  to  hold  to  get  a  sero’th  order  correction  of 


*J|  "  I  “J0  *  »  "_*  "  P  "  S  ’  •  ‘  >  *  »  •  i  •  »  '  i  *  »  "»  '  •  •  »  ■  •  »  *  •  »  •  •  •  •  ‘  *  a  '  •  -  t  k  •  . 

:  - :  ;  >:-v  - 

-•  •  . . . . . -  .wv  ,v.  -  . 


^  .• 


e  •  .  -  *  •  *  *  *  *  *  •  «  *  a  '  •  *  a 

v-v-*. 


.V  .*■  .V 


r% 


vyv. 


1+-+1 
I , 

p 

u 

Pi 


t 


V 


! 

/ 

> 

! 


fr. 


(5.20)  Proof  of  Theorem  (3.33);  Let  y  be  defined  by 

(5.21) 


—0  ax 

y(x,x)  -  4>  (x,x)e 


for  some  a  satisfying  (3.37).  Then  y  satisfies 

(5.22)  yxxT  -  -(<2p°(x,0)  +  C(x)  -  ct)yx)x 

(5.23)  y(0,x)  -  y(1,x)  -  0  . 
Multiplying  (5.22)  by  y  and  integrating  by  parts  gives 


(5.24) 


If 


2  3^  /  Vxdx  “  “  /  ( 2p(x,0)  +  C(x)-o)yxdx 


min|2p(x,0)  +  C(x)  -  a|  >  0 
x 


holds  then  the  L2  norm  of  y  stays  bounded  and  (3.25)  holds. 
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